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In order to eliminate disturbing effects of decoherence, encoding of quantum information in 
decoherence-free subspaces has been suggested. We analyze the benefits of this concept for a quan- 
tum register that is reahzed in a spin chain in contact with a common bosonic bath. Within a 
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master equation of Bloch-Redfield type. 
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I. INTRODUCTION 

The main obstacle in utilizing the remarkable 
computational power of quantum systemsii^i^ is 
the omnipresent and fundamental phenomenon of 
decoherence^!^. While this insight cast significant doubts 
about the idea of large-scale quantum computation^"^, 
it also initiated extensive research on decoherence in 
quantum information systems, and, beyond that, it led 
to the development of quantum error-correctin g^i^i ^'^ 
and -avoidingiiii^ii^ methods. The latter ones, on which 
we focus in the present work, make use of possible 
symmetries in the interaction of, say, a quantum reg- 
ister and its surrounding environment. The idea is to 
encode quantum information in those register states 
that are protected by symmetry against the decohering 
interaction. Inasmuch as the symmetry is satisfied these 
states span a decoherence-free subspace in the register's 
Hilbert space. 

Almost necessarily physical realizations of this con- 
cept will have to rely on symmetries that hold only to 
some approximation. Encoding in subspaces that respect 
these symmetries can then provide only partial protec- 
tion against decoherence, to an extent that will depend 
on the actual realization. The present work addresses 
this problem for the generic situation of a quantum reg- 
ister consisting of (effective) spin- 1/2 particles in contact 
with a common bosonic bath. 

In the main (Sees. Illl and IIIl| ). we describe this sys- 
tem by the dissipation-less spin-boson model of Palma 
et alJ^ which has been frequently used in similar 
contextsi^ii^ii^iii. For simplicity, we assume the spins 
to be arranged in a linear chain with inter-spin dis- 
tance a. Furthermore, we will use a three-dimensional 
bosonic bath with an ohmic coupling density of states. 
In the limit of vanishing distance a the model exhibits 
a highly symmetric spin-boson interaction, allowing the 
construction of decoherence-free subspaces. Specifically, 
we consider subspaces that correspond to encoded quan- 
tum registers in which logical qubits are encoded in lo- 
cally grouped physical qubits (spins) (cf. Sec. IIIB[) . 

The main task is to analyze the decoherence which will 
appear in these encoded quantum registers when the dis- 
tance a assumes finite values. In Sec. Illll we quantify the 



decoherence of encoded registers (with finite distance a) 
as well as of plain registers by means of the average regis- 
ter fidelity. Sec. llVl is devoted to the effect of dissipative 
spin-boson couplings. 

For a summary of our results we refer to the self- 
contained presentation in Sec. |Vl 



II. DECOHERENCE OF QUANTUM 
REGISTERS 

A. Physical quantum register 

A physical n-qubit register may consist of n (effective) 
spin- 1/2 particles located at sites ro, r2, . . . , r„_i of a one- 
dimensional lattice with lattice constant a and of finite 
length L = a(n— 1). A homogeneous (effective) magnetic 
field in z-direction may lead to a Zeeman energy splitting 
£. The corresponding register Hamiltonian is 

where Zi denotes the Pauli (Jz-operator for spin I. The 
register is supposed to be in contact with a thermal bath 
of three-dimensional bosons described by the Hamilto- 
nian 

k 

where and are creation and annihilation operators 
of bosonic modes with linear dispersion ~ c|k|. We 
assume a linear and local spin-boson interaction via Zib\^ 
and Z/6k operators. The corresponding coupling con- 
stant (7;^k will acquire a phase e^^'^\ reflecting the wave- 
like character of the bosonic modes. Apart from this 
phase the interactions may be isotropic and identical for 
each spin. Thus, g;_k = <?|k|e''''''' , resulting in an interac- 
tion Hamiltonian 

Hrb = Y.ZiB{vi), 
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where B(r) is the hcrmitian bosonic field operator 



with a distance and time dependent decoherence function 



i?(r)=^g|k|e-"^- 

k 



H.c. 



As customary, wc describe the strength of the spin-boson 
couphng by a spectral density 



auj e 



It is characterized by a cut-off frequency fl, a constant 
a of appropriate dimension, and a non-negative spectral 



parameter 



„18.19 



Since the spins interact with the bosons 



only via the energy conserving Zi operators, the model 
shows no dissipation but pure decoherence. This restric- 
tion makes the model analytically manageable F^iii^iiiii^ 
(Effects of additional dissipativc couplings will be dis- 
cussed in SeclIVn 



1. Decoherence of a physical quantum register 

The register may be used to store quantum informa- 
tion in form of a state po in which it is initially prepared. 
In general, po is subjected to a non-unitary dynamics 
originating from the system's own dynamics and its cou- 
pling to the bosonic bath. Assuming that the bath is 
initially in a thermal state ps, the total initial state is 
Qa ^ Pa ® Pb- During some time period t this state will 
evolve unitarily according to ihg = [Hfi + Hb + Hub, q] 
towards a final state gt- Its partial trace with respect 
to the bosonic modes yields the reduced density opera- 
tor pt = tr BQt, which describes the final register state. 
This procedure defines a quantum opcratio n^^i^*^ £ on the 
register by 

Po i-^ £{po) := Pt ■ 

The work of several author o^'"'^"'^i"'^^'"'^^i"'^''''"'^^i^"'^ established 
the operation £ to be of the form 

£{p)=UoAf{p), 

where U is a, purely unitary operation, and A/" is a non- 
unitary, completely positive map that can be written as 



■\p){p\p\l^){l^\ 



(1) 



Here, the double summation extends over all regis- 
ter eigenstates {v) which wc label in the usual way 
by n-bit sequences u G {0, 1}" = according to 
= {—ly'liy). The decoherence coefficients D^i, 

are 

n-l 

Df,^^ ^{pi-iyi){pm-iy,n)K{\ri-r,n\,t) , (2) 

lm=0 



Ki\rlt) 



4Re J dt' J dt" {B{r,t")B{O,0))T 



(3) 



Here {■ ■ ■)t denotes the thermal average over the bosonic 
system at temperature T, and B{r, t) is the bosonic field 
operator in interaction picture, i.e. 



B{r,t) 



k 



(4) 



The unitary part U of £ originates from the registcrs's 
own dynamics but also includes the Lamb-shift caused 
by the bosonic bath. In principle, this part of the evo- 
lution £ can be reversed and therefore is not of major 
concern. In contrast to that, the operation Af gives rise 
to decoherence and so seriously affects the register in an 
irreversibly manner. Clearly the main attention has to 
be paid to Af. Therefore, it will be in the focus of our 
investigation. 



2. Ohmic decoherence function 

In order to make Af more explicit, we have to deter- 
mine the decoherence function Eq. ([3]). For the sake of 
simplicity, we restrict ourself to an ohmic spectral func- 
tion J(cj) = awe""/^. Then, determining the correlator 
{B{r,t")B{O,0))T by standard methods and passing the 
continuum limit for the bosonic modes, the decoherence 
function ^ becomes 



K{r,t) = a / dw 



1 — cos ivt , / w \ sin ujr , , /o 

coth — e"'^/" . 

uj \2TJ ujr 

(5) 

(Henceforth we use units in which c = 1, ft = 1, and 
~ 1.) This integral can be better dealt with if we 
distinguish between the case of strictly vanishing distance 
r and the case of a finite distance that is large compared 
to the cut-off wavelength ^ which we assume to be 
the smallest scale in the problem. 
For vanishing distance we obtain 



K{0,t) = 2a In 



T{T/n) 



-Hi + t^n'). (6) 



T{T/n ~ itT) 
At small times i >C 1 /T this is in good approximation 



A^(0,t) -In(l + <2f72), 

whereas for large times t ^ 1 /T we have 

K{0,t) ~ QfTrTt + a In — — . 

Physically, K{0, t) determines the decoherence of a single 
spin, as it is seen by Eqs. ([2]) and ([T]) for n = 1, according 
to which 



(7) 



(8) 



Poi(t) 



-K{0,t] 



Poi(O) 
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FIG. 1: Decoherence functions K{0, t) and K{ro,t) as a func- 
tion of dimensionless time r = t/ro at temperature T = 5/ro. 
The cut-off energy is n ^ lO^T. The dotted and dashed 
curves are short- and long-time approximation, respectively. 



Asymptotically, the single-spin decoherence decays expo- 
nentially with a rate 7 = airT, by Eq. ([5]). 

For finite r ^ the oscillations of the spherical 
Bessel function sin(a;r)/a;r damp the integrand in Eq. ([5|) 
more effectively than the regular cut-off exp(— w/il). 
This allows us to take the limit fl ^ 00. The result- 
ing integral can be solved by contour integration, which 
finally leads us to 



Kir,t) 



m:T{t 
a 



1 



-) + 



inTr 



12T^r 

[fT{t + r)-fT{t-r)-2fT{r 



(9) 



for r < t, and to 



K{r,t) 



anT ■ 



2r 



4:TTTr 



[fTit + r) + fTir-t) 



(10) 
2. Mr)] , 



B. Encoded quantum register 

For vanishing lattice constant a the locations of all 
spins of the quantum register introduced in III Al fall onto 
a single point Fq. This implies a highly symmetric spin- 
boson interaction 



(13) 



that exactly annihilates all states with vanishing total 
spin-z component. As a consequence, any linear subspace 
C of the register's state space ff„ that is spanned by 
such states is not affected by the bosonic bath at all. It 
represents a decoherence-free subspac o^^'^^i^^ . 

At a finite lattice constant a the former symmetry is 
absent and consequently C ceases to be decoherence-free. 
However, by reasons of continuity the decoherence of 
states in C will be still much lower than for arbitrary 
states as long as the lattice constant a is not too large. 

Of course, the decoherence reduction at a finite lattice 
constant a will strongly vary for different choices of the 
subspace C. Here we will investigate subspaces that result 
from encoding (logical) qubits in local groups of physical 
qubits. This is supposed to be done in a regular manner 
such that the resulting structure forms a regular encoded 
quantum register. 

To be specific, let us consider a one-dimensional phys- 
ical 2n-qubit register i?2n whose 2n spins 5*0, ... , S'2n-i 
are grouped in n pairs of neighboring spins as sketched 
in Fig. [21 The spin pair S2iS2i+i has a four-dimensional 
Hilbert space of which the two orthonormal states 



|0>!:= 10)2,11; 



2'i+l, 



|1>2.|0)2.4 



(14) 



are annihilated by the spin-boson interaction if a = 0. 
In this case the subspace Ci spanned by states \0)j and 
is decoherence-free. We call Ci the state space of 
the encoded (logical) qubit Qi, and we further define an 
encoded n-qubit register to consist of the n encoded 
qubits Qo, . . . , Qn-i- Its Hilbert space = Co (Xi • • • (S) 
Cn-i is by construction decoherence-free with respect to 
^RB ■ ^il^ denote a state p as a state of the encoded 



for r > t. For convenience, we introduced a temper- register Rj, if the support of p lies entirely in 



ature dependent function frit) := Li2(e~ ), where 
Li2(x) = x-'/j^ is the dilogarithm of x. 

In the large temperature regime characterized by 
t,r, |t - r| > l/T, the /T(-)-terms in Eqs. ^ and ^ 
become exponentially suppressed, and so, additionally 
omitting an Q7r/24Tr term. 



K{r,t) 
K{r,t) 



aiiT (t 

2 



2r 



, forr<i (11) 
for r > i (12) 
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|1)° |0)° |0)° |0)° |1)° |1>° |0)° 

|o>i |o)i 

|1>= |0)= 



qubits 
1^* order encoded qubits 
2nd Qj-cier encoded qubits 



Fig. [T]shows K{0^t) and K{rQ^t) as a function of dimen- 
sionless time r = t/ro together with their approximations 
Eq.s (HD,®, CID, and for Q = lO^T. 



FIG. 2: Spin array representing qubits and encoded logical 
qubits of 1"* and 2"'' order. 

Clearly, pairing up remote spins instead of adjacent 
ones would lead to encoded qubits that will be more sensi- 
tive to an increasing lattice constant a. We will therefore 



4 



exclude this possibility from our considerations. Instead, 
one may speculate that the protection against decoher- 
ence improves if we iterate the pairing in order to built 
encoded qubits and registers of higher order (cf. Fig. [5]). 

More precisely, we recursively define encoded qubits of 
order x = 1, 2, . . . by their logical states 



mi 



io)r'ii)2T+i, 



ID 



ID 



+1 ' 



where 0*'' order states arc identified with plain spin states 
|0)? = |0)i and |1)° = \ of spin Si. An encoded n-qubit 
register of order x is then built from the encoded 
qubits of a 2n-qubit register i?2^ of order X ~ 1- ^ 
state of the encoded register i?^ is by definition a 
state whose support lies in H^. 



1. Decoherence of encoded quantum register 

How will the decohering operation M affect the en- 
coded registers which we have just introduced? First 
we observe that a state of an encoded register i?^ 
remains a state of under A/", simply because the spin- 
boson interaction Hub does not flip spins. Moreover, in 
App. El we show the operation M on an encoded register 
i?^ to be formally given again by Eqs. ^ and ([2|). What 
changes is the decoherence function K{\Yi—Ym\,t)^ which 
has to be replaced by an effective decoherence function 
Ky^_,{t), and the summation, which now extends over 



logical register states \p)^, \ v)^ given by 



Imo)^ . 



iMn-l/n-l 



Explicitly, for a state p^ of an encoded register we 
have 

M{p^)= e-^-|M)'^(M|Pxl^)'^(Ml, (15) 

with effective decoherence coefficients 

n-l 

Dlu = E - '^')('^'" - • (16) 

im=0 

The effective decoherence functions Kf{t) for x ^ 1 
recursively defined by 

Kf{t) ^ 2K^-'(t) - Kf-\{t) ~ Kl^-^\{t) , (17) 
with K^{t) = K{la,t). 

C. Discussion 

The formal analogy of Eqs. (US]), dm) and Eqs. 1^ © 
allows for a first comparison of the decoherence in en- 
coded and plain quantum registers by simply compar- 
ing the corresponding decoherence functions given by 
Eq. (HZI and Eq. ©. 
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FIG. 3: Decolierence functions K{0,t) (dashed) and Ko{t) 
(solid) as function of dimensionless time t = t/a at high (a) 
and low (b) temperatures with cut-off Q = lO'^/a, and T = 
10/a in (a) and T = 0.1/a in (b). 



We begin with the 1** order decoherence function 
KQ^t), which describes the effective decoherence of a sin- 
gle encoded qubit i?}. According to Eq. pT|) 

In the high temperature regime t,a^ 1/T we may use 
approximations ([5]), (fTTI) . and (fT^ to derive 



K^{t) = 2a In 



n 

2^ 



airTa : t > a 

aTTT{2t- ■^) : t<a 



The effective decoherence function (t) increases twice 
as fast with time as A'(0, t) for small times, but quickly 
saturates to a constant value at time t ~ a (cf. Fig. [3^). 
Qualitatively, this remains to be also true at lower tem- 
peratues. For a<Cl/T'<Ctwe can extract from relation 
(ISI) and the exact expression ([9]) an asymptotic value 

Kq{oo) ~ 2aln — 

that is reached again at t ~ a (cf. Fig. [31d). 

We conclude that for any finite distance a the co- 
herence e~^o'-*'' of an encoded qubit approaches a fi- 
nite asymptotic value at times t > a. In the long- 
time limit the coherence of the encoded qubit will there- 
fore largely exceed the exponentially decaying coherence 
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g-K(o,t) ^ (27rT/r2)"e-"''^* of a plain qubit. At short 
times, however, the encoded qubit performs worse than 
the plain qubit. The crossover time t(. can be easily de- 
termined to be tc — a in the high temperature regime 
T » 1/a, and ~ ^ In(27ra2f]r/e2) > a in the low 
temperature limit T <^ 1/a. Thus, whether it is bene- 
ficial to encode or not also depends on the time period 
over which the qubit is supposed to store information. 
Here, it is important to observe that with lowering the 
temperature T one eventually reaches the low tempera- 
ture regime where the crossover time tc increases with 
l/T. 

Do things further improve when one goes to higher- 
order encoded qubits? Interestingly, this is not the case, 
for the reason that the 1*** order qubits of an encoded 
register i?^ are already essentially decoupled (see below) , 
and pairing up these independent qubits to higher-order 
qubits would not further reduce their effective decoher- 
ence. The decoupling of the I''* order-qubits is seen from 
the effective decoherence functions Kl{t) for / > 1. By 
Eq. ini) we find 

Kl{t) = 2K{2la,t) - K{2la-a,t) - K{2la + a,t) . 

In the high temperature regime this predicts by Eq. (jll[) 
actually a vanishing Kl{t) for times t > 2la + a. More 
precisely, for / > 1 

Alternatively, we can also directly calculate the zero- 
distance decoherence functions -STq (t) according to rela- 
tion (jl7p . In the high temperature regime wc obtain in 
the long-time limit for x > 1 

K^{^) = 2^''Kl{^), 

which obviously strongly increases with x- Qualitatively 
similar behaviour is found also at low temperatures, 
where we used the exact result Eq. ([9]) to numerically 
determine Kq{oo). The results are plotted in Fig.|4] 

So far the discussion is restricted solely on a compari- 
son of (effective) decoherence functions. While this suf- 
fices to characterize the decoherence of (encoded) qubits, 
it does not necessarily provide full insight in the perfor- 
mance of an entire (encoded) quantum register under the 
noise operation A/". The next section addresses this point 
by a systematic investigation of average register fidelities. 

Another important point which deserves further inves- 
tigation is dissipation. The above analysis is based to 
a large extent on exact results that are available only 
for the dissipation-less spin-boson model. Therefore, it 
might be possible that essential features of encoded reg- 
isters - particularly the saturation of the decoherence - 
will not survive when dissipativc spin-boson couplings are 
taken into account. We will investigate this problem in 
SecHVl 
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FIG. 4: Asymptotic value of the decoherence functions 
Kq{oo) as function of dimensionless temperature Q = Ta 
for X = 1, 2, 3 and cut-off ^ = 5 ■ lOV^- 



III. AVERAGE REGISTER FIDELITIES 

A. Definition 

Let a physical or encoded quantum register i?„ with 
Hilbert space iJ„ be exposed to some noise operation M 
(for instance, the one studied above), under which an 
initial pure register state ip = l'0)('0l evolves to a final 
J\f{ijj). The channel fidelity of ijj with respect to JV, 

captures how well the state is preserved in this process^. 
A quantum register will be suitable for information stor- 
age under the noise Af if, in average, the channel fidelity 
F{ip,Af) for register states is large. A reasonable figure of 
merit is therefore given by the average fidelity of register 
Rn with respect to Af, 

Jd^ F(i^,J\f), N = y'dV' 1, 

where the integrals extend over all pure code states tp 
with state vectors Itjj) G Hn with respect to a unitary 
invariant measure. 

More precisely, we can express the average by an in- 
tegration over the group U(-ffn) of unitaries on Hn with 
the normalized Haar measure fi, 

F= J d^i{U) FiUi'oU\N) , (18) 

U(//„) 

where ■00 is any fixed pure register state. 

For large qubit number n the average register fidelity is 
known to agree with the entanglement fidelity of Af with 
respect to Hn , which plays a prominent role in quantum 
information theory. 
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B. General expressions 

1. Sum representation 

From now on we consider the noise operation A/" of the 
dissipation-lcss spin-boson model which wc introduced in 
Sccini Our aim is to derive an expression for the average 
fidehty for physical or encoded registers in terms of the 
decoherence coefficients ([T6|) . Using representation 

2"-l 
fiiy—O 

(we omit the order-index ^-nd identify with in- 
teger numbers 0, 1, 2, . . . , 71 — 1) and relation ([18]) with 
ipo = |0)(0| immediately obtain 

2"-l „ 
""^^ U(W,.) 

The integral can be calculated by standard methods^^, 
leading to 

U(H„) 

This outcome is easily understood once one recognizes 
the integration as an average over all 2"-dimensional 
complex unit vectors uq e Hn of which Uq^ and Uqi, are 
the /ith and vth component, respectively. For large di- 
mension 2" and [i ^ v the squared absolute values |?7o/jP 
and |t/oiyp are nearly independent gaussian variables X\ 
and X2 of mean (X\) = (X2) = 2~", by the normal- 
ization of Uq. For [i ^ V the integral over the product 
therefore amounts to the expectation value 
{X-yX-ij = (2~")^. For [I = V instead obtain the sec- 
ond moment {Xl) 2(Xi)2 2(2-")2. The exact calcu- 
lation reveals rather tiny corrections to these estimates of 
order 2~^", which we will neglect in the following. Then, 
by the last two equations we find 

^ = |rE-"'''"^(i + '5.'')- 

iiu—^ 

Since the sum over the extra diagonal terms er^'^"^p,v 
contributes at most 2~" to the fidelity we can omit these 
terms as well and thus are left with 

^ = ^E^-"-- (19) 

iiv—^ 

In general, the double summation over the exponentially 
large range 0, . . . , 2" — 1 makes a direct numerical or an- 
alytical evaluation of this expression difficult. However, 
progress can always be made if we proceed similar as 
in Ref.— and employ a Hubbard-Stratonovich transfor- 
mation. This will factorize the double sum into n trivial 
sums, at the expense of an n-dimensional integration over 
auxiliary continuous degrees of freedoms. 



2. Integral representation 
We rewrite the decoherence coefficients Eq. (|16p as 

where we introduced real, n-dimensional vectors v^^ with 
components 

and a real and positive n x n decoherence matrix K 
whose entries are determined by the (effective) decoher- 
ence functions Eq. (|17p . 

Kim = K\i_ra\{t) . 

Then, with Gauss's identity 

-v+„Kv^„ _ /•d"x xtK"^x + 2»vt„x 

where N = (7r"detK)'/^ the average fidelity Eq. (fT9|) 
becomes 

J A^N ^ 
The sum is readily determined to be 

2"-l ri 1 n 

such that the average fidehty becomes 

['^-^e-'^''-f[cos'x,. (20) 

1=1 

This relatively well-behaved integral representation of 
the average fidelity can serve as starting point for nu- 
merical or analytical calculations (cf. Sec. IIII C|l . Fur- 
thermore, in contrast to the sum representation (|19p . the 
integral (|20[) indicates how to obtain approximative ex- 
pressions. 

3. Weak coupling approximation 

For weak couplings a the inverse eigenvalues of K be- 
come large and hence the integrand sharply peaks at the 
global maximum at x = 0. In this case it is appropriate 
to expand the integrand as 

n 

g-xtK-x-Q^Qg2^^ = e-'"(K-Hl)x + 0(|x|^) 
1=1 

and to omit the 0(|x|^) corrections in the exponent. In- 
serting this in Eq. (j20[) wc arrive at a proper Gauss in- 
tegral which yields a surprisingly simple weak coupling 
approximation 

J^^c = det(l + K)-i/2 (21) 
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of the average fidelity. By consideration of a diagonal 
matrix K we estimate the relative error | F — Fy^c \/F oi 
order trK^, which one has to keep in mind when using 
this approximation. 

4-. Small deviations 

Finally, let us consider the practically relevant situa- 
tion where the average fidelity deviates only by a small 
amount e from unity. By the weak coupling approxima- 
tion we find 

F = 1 - £ = dct(l + K)-i/2(i ^ 0(trK2)) , 
and, taking the logarithms of both sides, 

ln(l -e) = -itrln(l + K) + ©(trK^) . 

When we expand the logarithms we observe that in lead- 
ing order e = ^trK -I- O(trK^), and hence, since the 
dccoherence matrix K has constant diagonal elements 

F =1 - ]^nKo{t) + 0(trK2) . (22) 

Small deviations of F from unity arc thus determined by 
the zero-distance decoherence function KQ(t), describing 
the decoherence of a single (encoded) qubit, and they 
grow linearly with the number n of qubits. 

C. Examples 

The following discussion of two illustrative examples 
will provide more insights in the average register fidelity. 
The results will be also useful in the subsequent compar- 
ison of plain and encoded quantum registers. 
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FIG. 5: Average fidelity of an 100-qubit register of indepen- 
dent qubits as a function of the single-qubit decoherence pa- 
rameter K — Ko (t) ■ The exact result (thick dashed curve) and 
the weak coupling approximation (solid curve) agree very well 
in the plotted regime. The upper dashed curve shows the ratio 
of exact and approximative fidelity. 



Not unexpected, the average fidelity of the 7i-qubit reg- 
ister is exactly the nth power of the average fidelity of a 
single qubit, Fi = (l-|-e~'*)/2. We note in passing that 
this result could have been derived also directly from the 
sum representation (|19p . 

We can also employ the weak coupling expression (|2ip . 
predicting 

F,,c = (1 + AJ)-"/^ . 

While this is not quite the exact result we in- 

deed observe good agreement for small couplings k <C 
^J\Jn. Notice that the weak coupling approximation 
particularly holds in the regime Xjn <^ n <^ yJTJn, 
where the average fidelity is already exponentially small. 

(cf. Fig. ED. 



1. Independent qubits 

The first example is a quantum register consisting of 
independent qubits, as it is reflected in vanishing deco- 
herence functions Ki{t) for / > 0. For instance, this is 
realized in a plain, physical register R'^ in the limit of a 
diverging lattice constant a, but also holds to good ap- 
proximation for an encoded register i?^ in the long-time 
limit (cf. discussion in Sec. Ill C|) . 

As a consequence of Ki{t) = for / > the decoher- 
ence matrix K of such a register is 

K = K 1„ , 

where k = KQ{t), and 1„ is the n x n unit matrix. Be- 
cause of this trivial matrix K the integral in (j20p nicely 
factorizcs into n one-dimensional Gaussian integrals. 



2. Symmetrically coupled qubits 

As the extreme opposite to the first, our second exam- 
ple is a register whose qubits are symmetrically coupled 
to the bosonic bath by an interaction Hamiltonian (jl3p . 
This is realized for a physical register in the limit of a 
vanishing lattice constant a, where all qubits are located 
at the same position. Consequently, here the decoher- 
ence matrix K becomes a uniform matrix with constant 
entries 

K,,„ = A'|,_,„| {t) = Ko{t) = K . (24) 

Up to a factor nn the matrix K describes the orthogo- 
nal projection on the diagonal d = (1, . . . , l)/y/n. K has 
therefore a non-degenerate eigenvalue nn with an eigen- 
vector d, and an {n — l)-fold degenerated eigenvalue 
with eigenspace d^. It follows that the integrand in 
PU)) has its entire weight on the diagonal d, as an effect 
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of which the ?i-dimensional integral collapses to a one- 
dimensional one. In this way the average fidelity results 
in 



1 



1 



2n 
1=0 



dx e 



-x-/nK ^Qg2n 



-K{n-iy 



The result can be better interpreted in the limit of large 
n ^ 1 and small k ^ 1 (independent of n) . When n is 
large, we can approximate the binomial factor 4~" (^") 
by a Gaussian, exp(— (72 — Z)^/n)/-\/7m, and further, when 
K is small, we are allowed to replace the sum by an inte- 
gral. This yields an average fidelity 



F = 



1 



(25) 



We notice that this expression also results from the 
weak coupling approximation Eq. (j2ip . since here 
det(l + K) = 1 -I- nK. 

The algebraical decay with n in Eq. ((25|) strongly con- 
trasts with the exponential decay of the average register 
fidelity Eq. p3|l observed for independent qubits. This 
marked difference must be attributed to the high degree 
of symmetry in the present case. In fact, the symmet- 
ric qubit-boson coupling (fT3|) entails that states with a 
small total spin-z component couple much less effectively 
to the bosonic bath as than they would do in the case of 
independent qubits. Apparently, for numbers of qubits 
and couplings with n 3> 1/k this results in a strongly 
enhanced averaged fidelity. 

Remarkably, a register with n 1/k does not bene- 
fit from these effects of symmetry. In this regime, the 
average fidelity for independent and symmetrically cou- 
pled qubits actually coincide (cf. Fig. [6]). We are lacking a 
simple explanation for that, however, since in this regime 
also 1 — F 1 we can refer to the general result Eq. p2|) 
for small deviations. According to this relation, here the 
fidelity is dominated solely by the zero-distance decoher- 
ence function Ko{t), and hence all details concerning the 
spatial structure of the register do not matter. 



D. Comparison of plain and encoded quantum 
register 

In this subsection we will compare a plain, physical 
register i?° with a 1** order encoded register i?^ by means 
of their respectice average register fidelities Fq and Fi. 
Thereby, we will make good use of the results for the two 
preceding examples. We will restrict the comparison to 
the high temperature regime a ^ where we can use 
the relatively simple expressions Eqs. (O, pT|) . and 
for the decoherence function. Furthermore, the time t 
will be assumed to be larger than Lq = {n ~ l)a for the 
plain register Rf and larger than Li — (2n — l)a for the 
encoded register i?^. 




0.00 0.02 0.04 0.06 0.08 O.IO 



FIG. 6: Average fidelities of an 100-qubit register for indepen- 
dent (solid curve) and symmetrically coupled (dashed curve) 
qubits as a function of the decoherence parameter At = Ko{t). 



1. Average fidelity of a plain quantum register 

We consider a physical n-qubit register (i?°) as defined 
in Sec. Ill Al In the high temperature regime and for 
t > Li, its time dependent decoherence matrix K(t) 
follows by Eqs. jS]) and (fTT|) to be given by 



Ki™(t) =7t I 



a\l 



2t 



7 = anT , (26) 



where we suppressed a logarithmic term 
a ln(ri/27rr) <C jt in the diagonal matrix elements. 
Since the non-trivial structure of K(i) does not allow 
for a simple evaluation of the exact formula Eq. (j20p . 
we immediately switch to a numerical evaluation of the 
weak coupling approximation Eq. (|2ip . 

Fo(t) =dct(l+K(t))-i/2_ 

The dashed curve in Fig. [7] shows the average fidelity 
Fo(t) of a linear qubit register with n ~ 125 qubits 
at a decoherence rate 7 = 10~*/a. The time domain 
130a < t < 1000a is chosen such that Eq. (gH]) and the 
weak coupling approximation is applicable. For compari- 
son. Fig. [7] also shows the average register fidelities of in- 
dependent and symmetrically coupled qubits with a time 
dependent parameter k = ^t, 



m = 



l + e- 



1 



(27) 



(cf. Eqs. §^ and ([25])). 

As expected, the average register fidelity Fo{t) lies be- 
tween Fi{t) and Fs{t). It might be more surprising that 
for large times < ^ Lq the fidelity Fg (t) is much closer to 
the fidelity Fs (t) of a register of symmetrically coupled 
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FIG. 7; Time dependency of the average fidelity for a linear 
quantum register of 125 qubits (dashed curve) in compari- 
son with the average fidelity of registers consisting of 125 in- 
dependent (bottom) and symmetrically (top) coupled qubits 
as a function of dimensionless time r = t/a. The rate is 
7 = 10-^ /a. 



qubits than to Fi{t), the fidelity corresponding to inde- 
pendent qubits. The reason for this resemblance of Folt) 
and Fs{t) is that for t ^ Lq the m-dependent term in 
Eq. ([26|) is a small correction to the leading term, mean- 
ing that K becomes close to the uniform decoherence 
matrix discussed in Sec. IIII C 21 with k = jt. 

At small times t < Lq, which are not covered by the 
decoherence matrix Eq. p6|) . relation Eq. ([22|) predicts 
the fidelity i^o(i) to agree with the converging fidelities 
for independent and symmetrically coupled qubits. 

2. Average fidelity of plain and encoded registers 

To begin with, we consider the average fidelity Fi of an 
encoded register at large temperatures T ^ 1/a. For 
times i > Li we find from relation Eq. (fTTl) with Eqs. ([8]), 
and a decoherence matrix 

Ku{t) - K^{t) = 7a, 7 = a^T, 
Ki^mit) = 0, 

where again we omitted a logarithmic term 
aln(ri/27rr) ^ 7a in the diagonal elements. Note 
that the time dependence has dropped out. This is 
precisely the decoherence matrix of a register with 
independent qubits which we analyzed in Sec. IIII C II 
Hence, by Eq. for times t > Li the average fidelity 
becomes a constant 

F,{^) = F,= ^—^ j • (28) 

For instance, for 7 = 10^^ /a and n ~ 125, which are the 
parameters of the fidehty Fo{t) plotted in Fig. ([7]), we 
obtain an -Fb(0 l^-rgely exceeding asymptotic value 

Fi(oo) ~ 1 - 0.0062. 



The saturation of the effective decoherence functions 
Kl{t) at times t > a entails the saturation of the av- 
eraged register fidehty Fi{t) to the asymptotic value 
i^i(oo). In contrast to that, the fidelity Fo{t) of the plain 
quantum register keeps decaying with increasing time. 
As a trivial consequence, the average fidelity of the en- 
coded register will always exceed the fidelity of a plain 
register if time t becomes sufficiently large. 

To resolve the time dependency of the average fideli- 
ties at shorter times we may use approximation Eq. ((^ . 
valid for small deviations I ~ F <^ 1, according to which 

Foit) = 1 - ^Koit) , F,it) = l- ^K'oit) , 

where Ko{t) and -ft'o(i) are decoherence functions for 
plain and encoded qubits, respectively. We therefore ex- 
pect encoding to be advantageous when Ko{t) > Kol{t). 
This is the case for times t larger than the crossover time 
tc which in Sec. Ill Cl was determined to het^ — a for large 
temperatures T ^ 1/a, and tc — ln(27ra^iir/e^) for 
low temperatures T ^ 1/a. 

We confirmed this by numerical calculation of the fi- 
delities -Fb(i) and Fi{t) withiir the weak coupling approx- 
imation Eq. (|2T|) . Fig. [8] shows the fidelity of encoded and 
plain registers in two temperature regimes (T ^ 1/a and 
T <C 1/a). The occuring crossover times arc in good 
agreement with the above discussed expectations. 

Finally, it might be instructive to analyze the require- 
ments on the single-spin decoherence rate 7 = anT 
which is needed in order to store quantum information 
within a given precision e during a given time period ts- 
(For simplicity, we restrict ourself to large temperatures 
T ^ 1/a.) Let us assume that ts scales with the number 
of qubits as 

ts ~ ton'^ , 

where q is some power greater or equal unity, and to ^ 
a is some microscopic time scale. Then, for the plain 
register the condition 

Foits) >l-e 

implies either by Eq. ((22|) or by Fo{t) w Fs{t) and 
Eq. ([27]) a rate 70 scaling as 

2s -(1+0) 
to 

On the other hand, for the encoded register the condition 
Fi > 1 — e is satisfied if, according to Eq. (^5]). the rate 
scales as 

2s -1 

71 n , 

a 

The advantage of using an encoded register over using a 
plain register is reflected in a possibly huge factor 

71 to „ 
— ^ — . 

70 a 
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FIG. 8: Average fidelity of encoded (solid line) and plain 
registers (dashed line) of 125 qubits for a large temperature 
T = 10/ a in the upper graph and for low temperatures 
aT = 1/5, 1/10, 1/20 in the lower graph as function of the 
dimensionless time r = '/a. The coupling is a = 10~"^/a. 
Notice that at low temperatures the fidelity of the encoded 
register is almost temperature independent. 



IV. DISSIPATIVE COUPLINGS 

The preceding sections have shown that the decoher- 
ence of an encoded qubit and also the fidelity of an en- 
coded register (of f** order) is essentially determined by 
the effective decoherence function 

In accordance with previous workiiiii we observed that 
for the dissipation-less model the decoherence coefficient 
KQ{t) saturates to a finite constant value Kq{oo) in the 
long time limit. It is important to find out whether 
this saturation holds for general physical interactions, or 
merely is a feature of the dissipation-less interaction. 

Here we cannot address this question in full generality. 
However, as a first step we will analyze a two-spin system 
that is dissipatively coupled to a bosonic environment. 
Lacking an exact analytical solution, we will treat this 
system by a quantum master-equation of Bloch-Redfield 
type (Sec. IIVX|) . Its viability in the present context is 
demonstrated for the dissipationless case, where the ex- 
act results can be reproduced (Sec lIVB]) . In Sec. lIVCl we 
will then consider the decay of a state that would remain 
invariant under the spin-boson interaction in the limit 



a — > 0. Contrary to what is observed in the dissipation- 
less case, here we find that for any finite a the state con- 
tinues to decay with a constant rate for large times. The 
asymptotic rate appears to be the given by the decay 
rate of a single spin multiplied with a universal factor 
2 - 2sin(£a)/(ea). 



A. Bloch-Redfield master equation 

We consider the general situation of a system S that is 
coupled to a bath B via an interaction Hj . The Lioville- 
von-Neumann equation for the density operator psB of 
the total system in the interaction picture is 

PS Bit) = -l[Hiit),psBit)]. 

It follows that the reduced density operator p (interaction 
picture) of the system obeys an equation of motioni^ 

m = - f dstrB[Hi{t),[Hi{s),psB{s)]] (29) 
Jo 

-iti-B[Hi{t),psBm- (30) 

We assume that initially ps b (0) is a product of an initial 
p(0) and a thermal bath state ps, and further take for 
granted that 



trs[ff/(i),p(0)®pB] =0, 



(31) 



which in many cases is satisfied, in particular in those 
to be analyzed below. The remaining term on the r.h.s. 
depends on the total state psBis) at times < s <t. In 
general, this does no allow to exactly determine p{0) from 
Eq. p9|) . One therefore frequently invokes the Born- 
Markov approximation by substituting 



Psb(s) ^ p{s) ® pb ^ p{t) ® Pb, 



(32) 



Obviously, this approximation is good in the limit of 
weak couplings. This results in the Bloch-Redfield mas- 
ter equation 



p(i) = 7^t(p(^)), 



(33) 



where TZt denotes the Rcdfield super-operator defined by 



p^Rt{p) = - dstvB[Hi{t),[Hi{s),p® pb]]. (34) 



Note that the Rcdfield operator is explicitly time- 
dependent and therefore the resulting dynamics does not 
exhibit a semigroup structure. In this sense, the Bloch- 
Redfield equation Eq. ((33|) is non-Markovian, notwith- 
standing the fact that the Born-Markov approximation 
has been used to derive it. In many cases it is justified 
to eliminate this "deficiency" by simply extending the 
domain of integration in Eq. from [0,i] to [— oo,t] 
(cf. Ref. fi^]). However, as it has been stressed by Doll 
et al Ji, when dealing with a spatially extended quantum 
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object this procedure would lead to noncausal behavior 
and thus to spurious results. 

In Ref. fi^] this problem has been circumvented by us- 
ing a causal master equation in which causality is explic- 
itly taken care of by step functions in the time domain 
that truncate acausal contributions. The resulting dy- 
namics has been shown to approximate quite well the 
known exact solution. Here, we will simply stay with the 
non-Markovian Bloch-Redfield equation as given by Eqs. 
dMl) and dsn). 

B. Dissipation-less two-spin system 

First, in order to demonstrate its viability, we use 
the Bloch-Rcdficld master equation to reanalyze the 
dissipation- less model of Section (III Ap for n = 2. The 
spin-boson Hamiltonian in the interaction picture is 

Hiit) = Zi®B(Yut), (35) 

1=0,1 

where Zi is the (time-independent) Pauli-z-opcrator on 
the l-th spin, and B{ri,t) as in Eq. ([4]). Condition ((3T1) 
is satisfied and the Redfield operator determines to be 

7^t(p) = ^d'^' - ^'"1' ^) (ZmpZi - ZiZ„,p) + h.c, 

m,i=0,l 

where 

t 

C7(|r|,t)- Jdt' {B{r,t')B{O,0))T . 



Presenting p{t) in the computational basis 

and again omitting imaginary parts, which would con- 
tribute only to the unitary pure U of the time evolution, 
the Bloch-Redfield master equation (|33ll predicts the co- 
efficients p^u{t) to obey independent differential equa- 
tions 

1 

Pliv{t) = ^ {p,l~Vl){p,rn~l'r,M^cC{c\m-l\,t)p^^{t). 
m,l=0 

After integration we observe that the master master- 
equation reproduces the exact result Eqs. ([2]) and ([3|) 
obtained in Sec. |lT] (for n = 2). This is more than one 
could have expected and we believe that the exactness 
must be ascribed to the fact that the present model is 
lacking dissipativc couplings. Wc do not expect that the 
Bloch-Redfield theory exactly describes the dynamics in 
the dissipative model. Nevertheless, the positive out- 
comes for the present dissipation-less model still encour- 
ages us to use the Bloch-Redfield master equation also 
for the dissipative model that we are investigating next. 



C. Dissipative two-spin system 

Let now the two spins interact with the bosonic bath 
via the dissipative Hamiltonian (interaction picture) 

i//(t) = ^ XKt)®S(ri,t). (36) 

1=0.1 

where B{r,t) is as in Eq. (|4]). Xi{t) can be conveniently 
written with operators 

ui = =Xi+ lYi 
di = a^!} =Xi- lYi 

as 

Xi{t)=e-'''di+e+'''ui. 

1, Master equation 

A straightforward calculation shows that in rotating 
wave approximation - which is valid as long as 3> 1 - 
the corresponding Redfield operator Eq. p4|) is given by 

T^t{p) = y^C-(|r; - rr,i\,t) {umpdi - diUmP} 

771,1 

+C+(|r; - 

771 

\,t){d 

771 

pui - Uid,nP} 

+h.c., (37) 
with energy dependent correlation functions 

C±(|r|,t)= / dse±*"(B(r,s)B(0,0)). 

The real part C'^ (r, t) of Cg (r, t) has a surprisingly simple 
structure in the limit of large Tt ^ 1. Here we find 
the zero-distance correlations to be given by the familiar 
expressionssi 

C^(0,i) = annie) 

C;(0,i) = a7r(n(e) + l), (38) 

where 71(e) = {e^^"^ — 1)^^. The finite-distance correla- 
tions follow to be (|r| = a) 

C^(a,t) ^ !!1^0(t_a)C;(O,t). (39) 

We notice that the finite-distance correlation deviates 
from the zero-distance correlation only by a temperature- 
independent factor that, of course, approaches unity for 
a ^ 0, but also for e ^ (when a < t), what is some- 
what unexpected. 

Anticipating the discussion given below, we remark 
that the deviation of the first factor from unity 

for any finite er reflects the decay of coherence in a sym- 
metric subspace in the long time limit. Note that for 
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e — > 0, which precisely corresponds to the transition to 
a dissipation-less model, the first factor remains unity 
also for finite a. Clearly, this corresponds to the satu- 
ration of the subspace fidelity observed in the previous 
section. The second factor, 9{t — a), which is responsible 
for causality, here as well as in the previous case leads 
to an initial drop of the subspace fidelity until the time 
t = a is reached. 



2. Asymptotic decay of subspace fidelity 

The anti-symmetric state 

l'/'o> = ^(|01)-|10)). 

is annihilated by Zq + Z\ and -f X\, and therefore 
remains invariant under both, the dissipationless inter- 
action Eq. ([35|) and the dissipativc interaction Eq. (pG]) . 
provided that a = 0. We are interested in the decay of 
the state |?/'o) under the dissipative spin-boson interac- 
tion Eq. ([36|) for finite distance a. To this end we consider 
the fidelity 

F{t) := (^oIpoWI'/'o), 

where /Oo(^) is the reduced spin state at time t (in inter- 
action picture) that evolved via the interaction Eq. (|36p 
with the bath from the initial state p(0) — |'i/'o)('0o|- 
Here we will approximately determine po (i) by the Bloch- 
Rcdficld master equation with the dissipative Redfield 
operator Eq. ((571) . 

To obtain a first impression of the dynamics we in- 
tegrated the master equation numerically. Characteris- 
tic outcomes for the fidelity F[t) are shown in Fig. [9l 
We chose energy e = 5T and a small overall coupling 
constant a = 0.01. Similar to the previously observed 
behavior, also here we see a relatively strong decay of 
the fidelity at times t < a. However, in contrast to the 
dissipation- less model, here we clearly see a decay of Fq (t) 
for large times with a rate that increases with distances 
ranging from aT = 0.1 to 0.3. 

At times t ^ a and for sufficiently small couplings the 
decay rate 71 = dF /dt of the fidelity F{t) is in good 
approximation determined by the expression 

71 =-(^o|7^^(|MV'o|)|^o>. 

Using expressions Eqs. (p8)) and (p9)) wc obtain after some 
algebra 



71 



2 1 



sin(ea) 



ea 



70 



where 70 = 2na{n{e) + 1/2) is the decay rate of a single 
spin. This is a simple and quite general result that - in 
the light of the discussions of the preceding sections - 
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FIG. 9: Fidelity of a symmetric subspace in a dissipative 
model as a function of dimensionless time t = tT with e — 5T 
for inter-spin distances a = 0.1/r(a), 0.2/r(fe) and 0.3/T{c). 
The larger the distance, the larger is the asymptotic decay 
rate a.t t > a. For comparison, the lower curve (d) shows the 
fidelity for the symmetric state (|01) -I- |10)). 



ro 




FIG. 10: The reduction 71/70 in the decay rate as a function 
of the dimensionless parameter p = ea. 



quantifies the benefits of using a symmetric subspaces in 
a dissipative system. It identifies 



P 



ea 
he 



as the relevant parameter that captures the achievable 
reduction of the decay rates of encoded qubits in com- 
parison to plain physical qubits. Fig. 1101 shows 71/70 as 
a function of p. 

A significant reduction requires p <C 1, corresponding 
to distances 

ch 

a < — . 

e 

For instance, for atomic qubits this implies that the dis- 
tance a between atoms should be small compared to the 
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wavelength of the light that is emitted in a bit-flip tran- 
sition. 



V. SUMMARY OF RESULTS 

In Sees . ITll and mil wc studied a quantum register phys- 
ically realized as a linear spin-chain that interacts with 
a three-dimensional bosonic bath via a dissipation-less 
spin-boson coupling with an ohmic coupling spectral den- 
sity. Within this framework we analyzed the benefits of 
using encoded qubits and encoded quantum registers in 
order to reduce effects of decoherence. 

In agreement with previous workii*^ we found that 
the coherence of a 1** order encoded qubit converges to 
a finite asymptotic value at times a/c, where a is the 
inter-spin distance and c is the velocity of sound or light. 
As a consequence, the coherence of the encoded qubit 
exceeds the one of a plain qubit for times larger than a 
crossover time tc- At high temperatures ksT 3> hc/a the 
crossover time is tc ~ a/c, whereas at low temperatures 
ksT <ti hc/a it increases roughly as tc — h/{kBT) with 
decreasing temperature T. Moreover, we observe the 1** 
order encoded qubits to be effectively decoupled from 
each other, meaning that higher-order encoding becomes 
counterproductive. On the other hand, this decoupling 
of encoded qubits should be advantageous for quantum 
error correction, which is known to be significantly ham- 
pered when the qubits are coupled via the bosonic batbi^. 
This aspect of using encoded qubits may deserve further 
investigation. 

In Sec, mil we derived a convenient integral representa- 
tion for the averaged register fidelity with respect to the 
noise of the dissipation-less spin-boson model. This al- 
lowed us to investigate the decoherence of an entire plain 
or encoded quantum register. Small deviations of the av- 
eraged register fidelity from unity are proportional to the 
number of (encoded) qubits and to the (effective) deco- 
herence function KQ{t) of a single (encoded) qubit. The 
improved performance of 1"* order encoded qubits there- 
fore carries over to an entire f** order encoded quantum 
register. This is confirmed by more detailed analytical 
and numerical results presented in Sec. IIIII 

Finally, in Sec. IIVI wc addressed the role of dissipation 
within a two-spin model. Its dynamics was determined 
by employing a master equation of Bloch-Redfield type. 
In the presence of dissipative spin-boson couplings the co- 
herence of an encoded qubit seems no longer to converge 
to a finite value. Instead, here we expect an exponential 



decay with an, however, reduced asymptotic effective de- 
cay rate 71 = 2(1 — sin(p)/p)7o, where the dimensionless 
parameter p is determined by the energy splitting e of 
the spins and the inter-spin distance a, p = ea/{hc). 
This work is supported by DFG grant No. KL2159. 
APPENDIX A: DECOHERENCE OF AN 
ENCODED QUANTUM REGISTER 

To begin with, let us consider a I*'* order encoded reg- 
ister R}^. Its Hilbert space is spanned by 2" state 
vectors 

Im)' = Imo)J . . . e H2„ , M e Z2 ■ 

Each p G Z2 corresponds one-to-one to a /i' G Z2" by 
demanding \p)^ = , which by definition ([T4|) means 

M2i = M» , M2i+i = 1 - M» ■ (Al) 

Let pi be a state of the encoded register, i.e. 
supp pi C . By Eq. ([T]) we find 

AA(pi)= e^"'"'lM)VlPi|^^)Vl, 
where ly' relates to ly as p' to p, and, by Eq. 

lm=0 

Making use of relation (|Aip and rearranging the sum we 
can rewrite D^'i,' in terms of p and v as an effective 
decoherence coefficient 

n-l 

= E - ''')(^™ - ^m)K\i_^^{t) , (A2) 

when effective decoherence functions Kl^{t) are defined 
as 

Klit) = 2KUt) - K\',,.^{t) - Rt^,{t) , 

with 

Kf{t)=Kila,t). 

It is straightforward to generalize this analysis to 
higher orders, which eventually leads us to relations (fTSj) . 
(fT6l). and (fT6l). 
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